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Abstract 

Background: The marine annelid Platynereis dumerilii is emerging as a powerful lophotrochozoan experimental 
model for evolutionary developmental biology (evo-devo) and neurobiology. Recent studies revealed the presence 
of conserved neuropeptidergic signaling in Platynereis, including vasotocin/neurophysin, myoinhibitory peptide and 
opioid peptidergic systems. Despite these advances, comprehensive peptidome resources have yet to be reported. 

Results: The present work describes the neuropeptidome of Platynereis. We established a large transcriptome 
resource, consisting of stage-specific next-generation sequencing datasets and 77,419 expressed sequence tags. 
Using this information and a combination of bioinformatic searches and mass spectrometry analyses, we increased 
the known proneuropeptide (pNP) complement of Platynereis to 98. Based on sequence homology to metazoan 
pNPs, Platynereis pNPs were grouped into ancient eumetazoan, bilaterian, protostome, lophotrochozoan, and annelid 
families, and pNPs only found in Platynereis. Compared to the planarian Schmidtea mediterranea, the only other 
lophotrochozoan with a large-scale pNP resource, Platynereis has a remarkably full complement of conserved 
pNPs, with 53 pNPs belonging to ancient eumetazoan or bilaterian families. Our comprehensive search strategy, 
combined with analyses of sequence conservation, also allowed us to define several novel lophotrochozoan and 
annelid pNP families. The stage-specific transcriptome datasets also allowed us to map changes in pNP expression 
throughout the Platynereis life cycle. 

Conclusion: The large repertoire of conserved pNPs in Platynereis highlights the usefulness of annelids in comparative 
neuroendocrinology. This work establishes a reference dataset for comparative peptidomics in lophotrochozoans and 
provides the basis for future studies of Platynereis peptidergic signaling. 
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Allatotropin, Neuroendocrinology, Proenkephalin 



Background 

Neuropeptides, including peptide transmitters and hor- 
mones, are a diverse group of signaling molecules involved 
in chemical communication among neurons and neuroen- 
docrine regulation. Neuropeptides represent by far the lar- 
gest group of neurotransmitters and neuromodulators [1] 
and are considered the oldest neuronal signaling mole- 
cules in metazoans [2]. Peptidergic signaling has deep evo- 
lutionary origins: essential enzymes for proneuropeptide 
(pNP) processing, maturation, and secretion have even 
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been identified in organisms without a nervous system, 
such as the sponge Amphimedon queenslandica [3] and 
the placozoan Trichoplax adhaerens [4]. 

pNPs are translated as inactive precursors typically 
consisting of an N-terminal signal peptide (SP) that di- 
rects the pNP to the secretory apparatus, and one or 
several peptide elements flanked by basic cleavage sites 
[5]. After pNP cleavage by neuronal prohormone con- 
vertases, the liberated peptides can be further modified. 
C-terminal alpha-amidation and N-terminal pyrogluta- 
mination are common forms of modification [6] and can 
affect peptide stability [7], peptide structure [8,9], and 
bioactivity [10,11]. 

Structurally and functionally important peptide elements 
often show sequence conservation among homologous 
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pNPs or between the mature peptides within a single pNP 
[12]. The strongest conservation is often restricted to a few 
key residues that confer bioactivity to the processed pep- 
tides [13,14]. pNPs often contain spacer sequences between 
the conserved peptide stretches. These spacers are less 
conserved and therefore considered to be biologically in- 
active [12]. 

Most neuropeptides signal via G-protein coupled re- 
ceptors (GPCR) [15]. It has become clear in recent years 
that GPCRs and their peptide ligands are stably associ- 
ated in evolution [16-18]. This co- evolutionary relation- 
ship of peptide-ligand pairs has been exploited to 
establish several conserved metazoan pNP families over 
large evolutionary distances [19,20]. 

The classical approach to identifying novel bioactive 
neuropeptides was biochemical purification from the 
species of interest, followed by functional analysis 
[21-24]. With increasing genomic and transcriptomic 
sampling, pNP identification has been accelerated by in 
silico sequence analyses based on homology to previ- 
ously described pNPs [25,26], or on the presence of se- 
quence features such as a SP, conserved C-terminal 
amidated motifs (e.g., RFa, "a" for amide) or cleavage 
sites [26-30]. 

Mass-spectrometry (MS) is also widely used as a 
powerful tool for the direct identification of bioactive 
peptides. This method relies on mapping the obtained 
peptide masses to a reference dataset (genome or tran- 
scriptome), and can be impeded by the presence of ex- 
tensive post-translational modifications [31-33]. A 
combination of genomics and MS approaches has re- 
vealed the complete neuropeptide repertoire of several 
species in many metazoan phyla [34-38] . 

Annelids represent a diverse and species-rich phylum 
and have long been used in neuroendocrinological and 
behavioral studies [39]. Comparative genomic ap- 
proaches [19,20,40] and other studies identified multiple 
annelid pNPs and neuropeptides, including RFa 
[13,41-47], FVRIa [48-50], excitatory peptide (EP) 
[51-53], egg-laying hormone (ELH) [54], myomodulin 
[55-57], RGWa [13], Lll or elevenin [28], vasopressin 
[39,58,59], gonadotropin releasing hormone (GnRH) 
[60,61], insulin-related peptides [62], neuropeptide Y 
(NPY) [63,64] and myoinhibitory peptide (MIP) [65]. 
Despite these advances, a complete picture of annelid 
neuropeptide diversity is still missing. 

Here we describe the neuropeptide complement of the 
marine polychaete annelid, Platynereis dumerilii, using a 
combination of transcriptomics, in silico pNP searches 
and MS -based peptide identification. Platynereis has re- 
cently been established as a powerful experimental or- 
ganism for comparative and experimental neurobiology 
[50,58,65,66]. Platynereis has a biphasic life cycle includ- 
ing free-swimming planktonic larval stages, followed by 



settlement and metamorphosis into the adult bottom- 
dwelling worm [67,68]. The Platynereis larval nervous 
system is highly peptidergic and several neuropeptides 
were shown to be involved in the regulation of larval be- 
haviors such as ciliary swimming and larval settlement 
[28,65]. Previous studies described 15 Platynereis pNPs 
that are expressed in specific neuronal populations. Pep- 
tides generated from these pNPs include various RFa/ 
RYa related peptides [28], vasotocin/neurophysin [58], 
FVRIa [50], RGWa [13], and MIP/allatostatin-B [65]. 

To complement this list, we used an integrative ap- 
proach and identified 98 Platynereis pNPs, most of them 
verified by MS analysis. Our pNP catalog represents the 
most extensive list of annelid peptides to date. This cata- 
log will provide a valuable resource for further studies of 
the peptidergic control of annelid behavior and physi- 
ology, and for the reconstruction of ancient metazoan 
peptide signaling systems and cell types [20]. 

Results 

Establishing the Platynereis transcriptome, predicted 
proteome and secretome datasets 

To identify novel pNPs, we performed deep sequencing 
of the Platynereis transcriptome using a combination of 
Sanger, Roche/454 and Illumina technologies. We se- 
quenced 77,419 expressed sequence tags (ESTs) from an 
arrayed, full-length normalized, mixed-stages cDNA li- 
brary [GenBank JZ391525 - JZ468943]. This library was 
further sequenced using the Roche/454 technology. We 
also obtained Illumina paired-end sequencing reads from 
13 Platynereis developmental stages including larvae, ju- 
veniles and adults. We assembled all acquired sequences 
into a reference transcriptome. The Platynereis tran- 
scriptome dataset contains 351,625 reads, with 87,686 of 
the contigs longer than 500 bp and 28,067 longer than 
1000 bp. The transcriptome was annotated using the 
Basic Local Alignment Search Tool (BLAST) against 
SwissProt and well-annotated bilaterian proteomes 
(Additional files 1, 2, 3 and 4). We also searched the 
transcriptome for open reading frames (ORFs) from 
which we derived a protein dataset. The predicted pro- 
tein dataset contained 51,767 sequences longer than 120 
amino acids (Additional file 5). To generate a Platynereis 
dataset of secreted proteins, the predicted protein data- 
set was analyzed for the presence of SPs [69]. We identi- 
fied 11,075 protein sequences with a SP. After the in 
silico removal of the SP, this secreted proteome database 
(Additional file 6), as well as the full predicted proteome 
database, were used for MS -based peptide identification. 

Identification of Platynereis pNPs 

In order to identify pNPs in the Platynereis transcrip- 
tome and predicted proteome datasets, we performed 
BLAST searches in these datasets using a large curated 
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set of metazoan pNP query sequences [19]. We also 
conducted pattern searches for repetitions of the motif 
x(3-10)-K[K/R]. The resulting sequences were examined 
for the presence of a SP, for cleavage sites, conserved 
peptide motifs, and other hallmarks of bioactive peptides 
and their processing (e.g. amidation signature C-terminal 
Gly, pyroglutamination signature N-terminal Gin, Cys- 
containing stretches, mono- or dibasic cleavage sites). 
These searches identified more than 80 Platynereis 
pNPs, including those previously described. 

To complement the bioinformatics screen, and to find 
evidence for the presence of the predicted active pep- 
tides, we performed liquid chromatography - tandem 
mass spectrometry (LC-MS) on peptide extracts from 
various Platynereis larval and juvenile stages. We 
mapped the obtained MS hits to the Platynereis datasets. 
Using the MS analyses we discovered 15 further pNPs 
that were missed during our bioinformatic searches, ex- 
tending the Platynereis repertoire to a total of 98 pNPs 
(Figure 1 and Additional files 7, 8, 9). We also analyzed 
the Platynereis neuropeptides for some of the most com- 
mon peptide modifications and detected several modi- 
fied peptides (Figure 1 and Additional files 8 and 9). 

All full-length Platynereis pNPs have a SP and at least 
one potential basic cleavage site, and they lack non- 
neuropeptide protein domains as defined in the PFAM 
database. Besides, all 98 pNPs had to meet at least one 
of several criteria to be considered as bona fide prohor- 
mones. These include homology to known metazoan 
prohormones (e.g., NPY, AKH, 7B2 [70,71]), a confirmed 
expression in Platynereis neurons (e.g., FVMa, SPY [28]), 
MS evidence (e.g., NGEW, GYa), conservation across 
lophotrochozoans or annelids (e.g., CCWa, QSGa) and 
other structural and functional hallmarks of a prohor- 
mone, such as a repetitive structure or peptide modifi- 
cations (e.g., SPRa, QRIa). 

Overview and classification of Platynereis pNP diversity 

For most Platynereis pNPs we obtained the full-length 
open reading frame sequence (including a SP) and for 51 
pNPs we have an available EST clone or a PCR-cloned 
cDNA. We annotated all pNP sequences with various se- 
quence features, including the presence of a SP, cysteine- 
rich stretches (potentially also involved in processing 
[72]), prohormone-convertase cleavage- sites, modified 
(amidated or pyroglutaminated) or non-modified active 
peptides and the presence of a MS -hit. For all Platynereis 
pNPs containing repetitive peptides we also generated 
sequence logos (Figure 1 and Additional file 8). We clas- 
sified the Platynereis pNP families according to their 
phylogenetic distribution, distinguishing ancient eume- 
tazoan, ancient bilaterian, and ancient protostome fam- 
ilies, as well as annelid/lophotrochozoa specific pNPs, 



and pNPs with no currently recognizable homologs 
outside Platynereis (Figure 1). 

For the phylogenetic classification, we performed a 
sequence-similarity-based (BLAST p-values) clustering 
approach [73]. We used a curated dataset of 6,225 pNPs 
from 10 phyla [19] combined with all Platynereis pNPs 
and their lophotrochozoan homologs collected from 
EST databases. We clustered this dataset using PSI- 
BLAST with 3 iterations. All Platynereis pNP sequences 
with no similarity to known metazoan pNPs were re- 
moved from the map. Metazoan pNPs that did not con- 
nect to any Platynereis pNP were also removed (e.g., 
parathyroid hormone, growth hormone) (Figure 2). 
Many of the repetitive pNPs formed a strongly con- 
nected cluster at the center of the map. These sequences 
were reanalyzed with clustering using non-iterative 
BLAST (Figure 3). The resulting maps were used to ob- 
tain an initial overview, which indicated a relationship of 
several Platynereis pNPs to known pNP families. The 
conservation of some Platynereis pNPs is limited to 
small stretches (the mature peptides) in the precursor 
and is difficult to identify using BLAST-based clustering 
(e.g., diuretic hormone 44 (DH44)). For this reason, we 
also performed motif searches and multiple-sequence 
alignments, which reinforced the family assignments, 
obtained by our clustering approach (Additional file 10). 

Careful inspection of unassigned Platynereis pNPs for 
the presence of short conserved motifs using the motif 
discovery tool MEME [74] and known peptide profiles 
[20] led to the assignment of further Platynereis pNPs to 
known families (e.g., DH44, SIFa, pyrokinin/small cardio 
active peptide (sCAP); Additional file 10). We also dis- 
covered several Platynereis pNPs that belong to newly 
defined lophotrochozoan and annelid pNP families 
(Figure 1 and Additional file 10). 

Platynereis pNPs belonging to ancient eumetazoan 
families 

Several Platynereis pNPs belong to ancient eumetazoan 
families [19] such as insulin-related peptide pNPs (IRPs) 
and the glycoprotein hormones bursicon A-l, A-2 and B 
(Figure 1). 

Platynereis also has three pNPs that belong to the eume- 
tazoan R[F/Y]a family, FMRFa, RYa, and YFa (Figure 1 
and Figure 3 and Additional file 10). The FMRFa pNP rep- 
resents the ortholog of other protostome FMRFa pNPs. 
RYa has orthologs in other lophotrochozoans, sharing 
a Pro-rich C-terminal peptide (Additional file 10). The 
Platynereis YFa pNP is also part of the FMRFa cluster 
(Figure 3), however, direct orthologs of YFa pNPs could 
not be identified outside Platynereis. 

A member of the eumetazoan Wamide/MIP/GWamide 
family, the Platynereis MIP/allatostatin-B pNP, gives 
rise to peptides involved in the regulation of larval 



Conzelmann et al. BMC Genomics 2013, 14:906 
http://www.biomedcentral.com/1471 -21 64/1 4/906 



Page 4 of 1 5 



fii*i 



R[F,Y]amide- 
Wamide/M I Pi'G Wamide - 



insulin-related — 



prokinelicin/astakm« 
glycoprolem bela 
vasopressin i' oxytocin 



RYamids {loptKHrodt > 
YFamda ^Pfaiynef&s} 
MIF alla"ostaLin-B 
Wl {.inmHiciJ 
IRP-1 



luqin 



allalotcoptn 
sCAP 
sutfakinin-1 
sulfakinin-2 
CCAP 

allalostalin-A-t 
gaianm/al laics tat m -A — | allatostalm -A-2 
allalostalmA-S. 
alleles tat*vC 
SIFam»de-1 
SlFamide-2 
L11| 

EP(j 



EomaloGtatin/allatostalin-C 
NPFF'GnlH/SIFamide-| 

NMB/CCHamide' Lr EP. 



IP 

| ■ 

C a 

HI 

Hi 



pdl 

laucokinin 



1 



myosuppressin/myomod. myomoduNn 



whithin/proctolin— | 



1 

wtirtmin-2 





Follcln-I 2*"£» „ 




1 FVamiCte/EFLGamicte 


I 


NKY I 




NKY-2 




FVRIamide 


U 
I 


GNXQN 


c 

B 


LXRX 


q 


CLCCY 


lopr 


CCWantta 


QSGamide 


CCRFannidft 
HFAarmda 


!5 


FVMafnkfe 


DLuiidt 


"3 


SL RFarriCe 


c 


QUEFtAS 


c 
(Q 


MNC 
LEQ 




WL □ 




SPY 




HK3A 




NGEW 




AYNPV 




YTL 




PPLPamida 




PY 


49 


Qpeplm 


i 


THDarrada 




YYQamide 


i 

1 


FGamide 


LPWamide 


GYarrvKte 


5 


SLL-1 
SLL-2 
SHM 
Rf4T 

SPRarrtHie 
QRIarridu 
QDNamtda 
FCNarnide 
CCAamide 



2 

IRP-3 
IRP-4 
IRP-5 

bursioon-A-1 
bufsioon'A-2 
bfunioon-B 
prokineticin 
giycoprotem-bela 
vasotocin -no urtphvsm 
NPY1 
NPY-2 
NPY-3 

RGWamide {toptxurDcn.) 

AKH-1 

GnRHJcorazonin'AKhURPCH— | AKH-2 
GnRH-1 

QnRH-2 
7&2 

proenkephahn 
CRF^rocortin^DH44-| DH^ 1 

achalin 
luqtn/RYamide 
tachykinin _| 

calciloniin/DH31— j ^ 

I pedal peptide 1 
orcokinin.'pcdal— I pedal' peptide 2 



orexin/allatol ropin 
nmed-U P BAIWpyrokinin 

eholacystokininfeulfakinin — | 

NPS/CCAP^conoCAP 



Figure 1 Neuropeptide precursor complement of Platynereis 

dumerilii. Neuropeptide precursors (pNPs) were classified based on 
their phylogenetic distribution into eumetazoan, bilaterian, 
protostome (prot), lophotrochozoan, annelid and P/afynere/'s-specific. 
Previously established metazoan pNP families are indicated in bold 
red [19]. For each pNP we indicated whether a cDNA clone, a full-length 
ORF sequence, and MS evidence are available. The presence of 
amidated, Cys-rich, or pyroglutaminated peptides, or a repetitive 
pNP structure, are also indicated. If a pNP family has multiple 
Platynereis members, we also indicate the likely origin of this paralog 
(e.g. {protostome}). 



settlement [65]. Platynereis has another related pNP, 
yielding non-amidated W[I/L] peptides (Additional file 
8 and Additional file 10). The two Trp residues that 
WI peptides share with MIPs (x-W-x 6 _ 7 -W-[G/I/L]), 
and the position of the WI pNP in the Wamide/MIP/ 
GWamide cluster (Figure 3) supports a close relation- 
ship between MIPs and WI pNPs. We only identified 
an orthologous WI pNP in the distantly related annelid 
Capitella teleta, suggesting that WI pNPs are annelid 
divergences of the MIP family. 

Platynereis also has a pNP with a cysteine-rich proki- 
neticin/colipase domain [75]. Directly after the SP, the 
Platynereis prokineticin pNP contains amidated LFVa 
peptides. A similar peptide could also be identified in 
the C. teleta prokineticin pNP (Additional file 10). 

Platynereis pNPs belonging to ancient bilaterian families 

Recent comparative genomic analyses defined more than 
25 pNP families that are ancestral to bilaterians [19,20]. 
We found Platynereis representatives of most of these 
families that could readily be classified as one-to-one 
or many-to-one orthologs. These include glycoprotein 
hormone beta, vasotocin-neurophysin, NPY, RGWa, 
adipokinetic hormone (AKH), GnRH, 7B2, proenke- 
phalin, DH44, achatin, luqin, tachykinin, calcitonin, di- 
uretic hormone 31 (DH31), pedal peptide, allatotropin, 
sCAP, sulfakinin, crustacean cardio-active peptide (CCAP), 
allatostatin-A, allatostatin-C, SIFamide, Lll, EP, pigment 
dispersing factor (pdf), and leucokinin (Figures 1, 2, 3 and 
Additional file 10). Two of these families, calcitonin/DH31 
and DH44, show interesting gene or peptide repeat dupli- 
cations in their evolutionary histories and will be discussed 
in detail below. 

Platynereis has two pNPs belonging to the calcitonin/ 
DH31 family, calcitonin and DH31. The calcitonin pNP 
shows high sequence similarity to the Cys-containing 
calcitonin peptides from vertebrates (Figure 4). In C. 
teleta, a second calcitonin-like pNP has recently been 
identified that is lacking the Cys residues [20]. We also 
found an orthologous Platynereis sequence. These annelid 
pNPs are more similar to insect DH31 pNPs, also lacking 
the Cys residues, as shown by clustering (Figure 2) and 



Conzelmann et al. BMC Genomics 2013, 14:906 
http://www.biomedcentral.com/1471 -21 64/1 4/906 



Page 5 of 15 



Platynereis * 

Vertebrata • 
Tunicata ▼ 
Cephalochordata ▼ 
Ambulacraria ▲ 
Annelida ▼ 
Mollusca a 
Platyhelminthes + 
Arthropoda • 
Nernatoda 
Cnidaria • 
Placozoa • 



83^ 



p-value 
1e-150i 

1e-10 



PBAr ^uqin 



IPY T rf / 49 T 



AT orexm 
T39 



glycoprotein- beta 

154 

glycoproteip-alpha 




pacapA/IP 



DH31 

CCAP-^43 

galanin ^ 



prohormone^ 



proctolin/ (Jrokineticin/ 
astakine 



T whitnin 



AstC 



16^ 

vasotocin/ 
neurophysin 



insulin- related 



Figure 2 PSI-BLAST cluster map of Platynereis pNPs and their metazoan pNP homologs. Nodes are colored based on taxonomy. Edges 
correspond to BLAST connections of P value < 1 e-1 0. Platynereis pNPs are highlighted as red stars. The identifier of Platynereis pNPs is provided in 
Figure 1 (CLANS-ID). 



sequence alignments (Figure 4). These results suggest that 
the ancestral bilaterian calcitonin peptide contained two 
Cys residues. The gene of this pNP duplicated in stem 
protostomes to give rise to DH31, followed by the loss of 
the Cys residues. Calcitonin has been retained in mollusks 
and annelids, but lost from the arthropod lineage. DH31 
was lost from mollusks, and preserved in annelids and ar- 
thropods (Figure 4C). 

We also identified two members of the corticotropin re- 
leasing factor/DH44 pNP family. DH44 has been shown 
to be related to mollusk egg-laying hormone (ELH) [20]. 
We identified two pNPs in Platynereis that showed high 
sequence similarity on the level of mature peptides to 
mollusk ELH and also to insect DH44, with several highly 
conserved amino acid positions (Figure 5). The Platynereis 
DH44 pNPs are highly repetitive (13 and 16 copies) com- 
pared to their mollusk or insect counterparts that have 
only one peptide copy per pNP (Figure 5). 



Platynereis pNPs belonging to ancient protostome families 

Several members of ancient protostome pNP families 
are present in Platynereis, including myomodulin, 
prohormone-3, and whitnin (Figure 1). 

The Platynereis myomodulin pNP has two alternative 
transcripts, both yielding several peptides with an LRMa 
C-terminus, characteristic of myomodulins (Additional file 
8 and Additional file 10). Comparison of the Platynereis 
myomodulin pNP with the C. teleta sequence revealed 
other highly conserved peptide-stretches, flanked by 
monobasic cleavage sites, which may potentially yield 
non-amidated peptides (PRXGX, Additional file 8 and 
Additional file 10). 

Lophotrochozoan-specific pNP families 

We identified various pNPs belonging to lophotrochozoan 
pNP families including the fulicin related pNPs FLa and 
FVa/EFLGa, NKY, FVRIa, GNXQN, LXRX, CLCCY, 



Conzelmann et al. BMC Genomics 2013, 14:906 
http://www.biomedcentral.com/1471 -21 64/1 4/906 



Page 6 of 15 



Platynereis * 

Vertebrata • 
Tunicata t 
Cephalochordata ▼ 
Ambulacraria a 
Annelida t 
Mollusca a 
Platyhelminthes + 

prepronoeiceptin 

Nematode 
Cnidaria • 

Placozoa ® ^ff* 



+ myosuppressin 



v. 



sNPF 

•Hi* i 



proenkephalin 



27 



POMC 



GnRH 



p- value 
le-2801 
1e-5 



55 



FMRFa ? 

« npf : 



62 



• «» • 



17 



,9 + 70 
fulicrn T 

TRH 4 * 



Ast-A 



44-46 



NPY 



3^ 



orcokmin 



MIP r pedai^i 
53 * jiedaW 



Nplpl 



leucokinin , , , . 

tachykir 



achatin 



.; akh « * » 



tachykinin 
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CCWa, QSGa, CCRFa, and HFAa (Figure 1). Many of 
these could only be identified in annelids, mollusks, 
and sometimes also platyhelminths. We refer to these 
as lophotrochozoan families for simplicity. However, 
further sampling in other lophotrochozoan phyla will 
be needed to clarify their history. 



Annelid FLa and FVa pNPs are related to mollusk fuli- 
cins, forming a lophotrochozoan pNP family. The Platy- 
nereis FVa mRNA has an unusual structure. Following the 
stop codon at the end of the coding region containing the 
FVa peptides, a second putative coding region is present, 
potentially yielding fulicin-like EFLGa peptides with an 
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Figure 4 Evolution of the calcitonin/DH31 pNP family. (A) Truncated multiple sequence alignment of Platynereis DH31 with arthropod DH31 
pNPs, all lacking cysteine residues. (B) Truncated multiple sequence alignment of Platynereis calcitonin with other lophotrochozoan and 
deuterostome calcitonin pNPs, all containing two conserved Cys residues. (C) Scenario of gene duplications and losses to account for the present 
distribution of calcitonin and DH31 pNPs in bilaterians. 
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Figure 5 Platynereis has two highly repetitive DH44 pNPs. (A) Precursor schemes of the Drosophilo DH44, Aplysio ELH, and Platynereis DH44-1 
and DH44-2 pNPs with the signal peptide (blue), cleavage sites (red) and predicted amidated peptides (yellow) indicated. (B) Multiple sequence 
alignment of protostome ELH and DH44 peptides. (C) Sequence logo generated from 16 different Platynereis DH44 peptides. 



extra Gly residue (Figure 6A). A similar structure has been 
described for the FVa/EFLGa pNP in C teleta [40]. The 
region encoding the EFLGa peptides lacks a start Met and 
a SP, therefore it is not clear whether this region could 
be translated to yield mature peptides. The analysis of 



Illumina reads did not reveal any alternative transcripts 
with in-frame EFLGa peptides following a start site and a 
SP. The conservation of the EFLGa peptide-stretches 
across annelids (Additional file 10) suggests that these 
peptides are functional. Since MS evidence for EFLGa 
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Figure 6 Neuron-specific translational stop-codon read-through of the Platynereis FVa/EFLGa pNP gives rise to two types of 
neuropeptides. (A) Precursor scheme of the Platynereis FVa/EFLGa pNP with the signal peptide (blue), cleavage sites (red), stop codons (grey) 
and predicted FVa and EFLGa peptides (yellow) indicated. (B) Anterior view of 48 hpf Platynereis larva stained with antibodies against EFLGa (left) 
and FVa (right). (C) Anterior view of 72 hpf Platynereis larva stained with an antibody against EFLGa (left), an antibody against EFLGa blocked with 
synthetic EFLGa peptide (middle) and an FVa antibody (right). (D) Dorsal view of 72 hpf Platynereis larvae stained for EFLGa (left) and a close-up 
of the EFLGa expressing neuron pair (left). (E) Dorsal view of 72 hpf Platynereis larva stained for FVa (left) and a close-up scan of a FVa expressing 
neuron pair (left). All samples were counterstained with a mouse acetylated-tubulin antibody (white). Green asterisks label neurons where FVa and 
EFLGa peptides co-occur. Scale bars: 50 urn. Abbreviations: pt: prototroch; mt: metatroch; patl: first paratroch; pat2: second paratroch; tt: telotroch. 
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peptides was missing, we raised and affinity purified a spe- 
cific antibody against EFLGa and performed immunos- 
tainings on Platy nereis larvae. The EFLGa antibody 
labeled two neurons in the dorsal episphere of 48 h and 
two pairs of neurons in 72 h larva (Figure 6B and C), in 
the region where a subset of the precursor-expressing cells 
are found [28]. When we pre-incubated the antibody with 
5 mM EFLGa peptide, we no longer observed cellular 
staining (Figure 6C). In agreement with previous studies 
[13,28], a cross-species reactive FVa antibody [13] labeled 
more neurons in the larval episphere at 48 h and 72 h 
(Figure 6B and C). EFLGa-positive neurons were a subset 
of these FVa expressing neurons as confirmed by their pos- 
ition, cell shape, and the presence of characteristic sensory 
dendrites abutting the cell bodies. (Figure 6D and E). These 
results suggest a translational stop-codon read-through 
mechanism [76] to yield bioactive EFLGa peptides, occur- 
ring in a subset of the precursor-expressing cells. 

Annelid- and Platynereis-spec\f\c pNP families 

Some Platynereis pNPs have recognizable orthologs only 
in annelids, including the FVMa, DLa, SLRFa, QERAS, 
MNC, and LEQ pNPs (Figure 1 and Additional file 10). 

For several pNPs we could not find orthologs outside 
Platynereis. Although we could not rely on pNP conser- 
vation as an indication of biological activity, we are 
confident that these pNPs produce bioactive peptides. In 
the case of WLD and SPY an effect on ciliary beating 
has previously been established [28]. Other pNPs harbor 
two or more copies of similar peptides, often confirmed 
by MS evidence, including HIGA, NGEW, AYNPY, YTL, 
PPLPa, PY, Qpeptin, THDa, FGa and SPRa (Additional 
file 8 and Additional file 9). Many pNPs contain peptides 
with hallmarks of biological activity including amidation 
(SPY, HIGA, AYNPY, PPLPa, THDa, YYQa, FGa, LPWa, 
GYa, SPRa, QRIa, QDNa, PCNa and CCAa), pyrogluta- 
mination (YTL, PPLPa, SLL-1, QRIa, QDNa) or various 
Cys residues that potentially form disulfide bridges (WLD, 
AYNPY, FGa, SLL-1, SLL-2, SHM, RNT, QRIa, QDNa, 
PCNa and CCAa) (Figure 1 and Additional file 8). Further 
sampling in annelids could lead to the identification of 
orthologous pNPs in closely related species. 

Stage-specific profile of Platynereis pNP expression 

The Platynereis transcriptome dataset was acquired from 
various larval, juvenile and adult stages. To profile pNP 
expression, we mapped the obtained stage-specific Illumina 
paired-end reads to all pNP transcripts. The total number 
of pNP reads increases through larval development, peaks 
in 15-day-old juveniles, followed by a drop in adult stages. 
There are also large differences in pNP expression levels 
between sexually mature males and females (Figure 7). 

We also performed a Pearson-clustering of pNPs 
based on their normalized stage-specific expression 



values. The Platynereis pNPs formed distinct clusters 
with expression peaking in different life cycle stages. 
We also found variable expression between 1 -month 
pre- and 1-month post-cephalic-metamorphosis, and 
sexually mature male and female samples (Figure 8 and 
Additional file 11). For example, FMRFa, allatostatin-C 
and allatotropin were highly expressed in males, and 
lowly expressed in females. 

Discussion 

An integrative approach to obtain the Platynereis pNP 
repertoire 

To obtain a broad complement of Platynereis pNPs, we 
used an integrative approach combining transcriptomics, 
peptidomics, and bioinformatics. Such an integrative 
methodology was indispensable to extend the Platynereis 
pNP repertoire to 98 sequences. 

In particular, MS analysis was often necessary to reveal 
pNPs with more restricted phyletic distribution, lower 
sequence conservation, or non-repetitive precursor struc- 
ture. For example, in the GYa and SLL pNPs only one short 
potential bioactive peptide occurs in the precursor sequence. 
Moreover, MS evidences combined with the conservation of 
the peptide in orthologous annelid and mollusk pNPs 
allowed us to conclude that these sequences represent bona 
fide pNPs with likely biological activity (e.g. GNXQN, LEQ). 

Expression profiling of Platynereis pNPs 

The transcriptional profiling of pNP expression through- 
out the life cycle revealed that most of the pNPs are 
expressed at multiple stages. This indicates that the 
Platynereis nervous system is strongly peptidergic 
throughout the life cycle. The broad temporal expression 
of most pNPs also suggests that several of the Platynereis 
neuropeptides have pleiotropic functions. For example, we 
have recently shown that 12 different neuropeptides affect 
ciliary swimming and larval settlement in early larval 
stages [28,65]. Since larvae settle after a brief planktonic 
stage (approximately 6 days) and locomotor cilia disappear 
after cephalic metamorphosis (approximately 1 month), 
the later expression of these peptides implies functions 
unrelated to ciliary swimming. The other specific differ- 
ences we uncovered between different stages may be re- 
lated to changing physiology (e.g., initiation of feeding). 
The differences between males and females indicate the 
presence of sex-specific neuroendocrine regulation po- 
tentially related to sex-specific physiology or behavior 
(e.g., in pheromone production or during the nuptial 
dance of Platynereis [68]). 

Platynereis, a powerful marine invertebrate model for 
studying the function of ancient peptidergic signaling 

Nervous systems evolved in a marine environment. Con- 
sequently, comparative studies of neuropeptide signaling 
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Figure 7 Stage-specific transcriptional profile of pNP reads. The sum of the total normalized reads per million (rpm) for the 98 Platynereis 
pNPs are shown for the indicated stage. The read counts for each pNP are normalized to the total amount of reads for each stage. 



in bilaterians would benefit from marine models with a 
rich set of conserved pNPs. Among the lophotrochozo- 
ans, a predominantly marine super-phylum, there are 
many emerging model species [77] . However, few studies 
used a comprehensive approach to identify pNPs in 
lophotrochozoans [20,25,40]. A recent study in the fresh- 
water planarian Schmidtea meditterranea established the 
first broad pNP resource in a lophotrochozoan [34]. Com- 
pared to S. meditterranea, Platynereis has a richer reper- 
toire of conserved pNPs. For example, S. meditterranea 
has no recognizable orthologs of achatin, luqin, AKH, alla- 
totropin, CCAP, pdf, leucokinin, tachykinin, proenkepha- 
lin and whitnin [34,78], all present in Platynereis. Previous 
studies also highlight the potential of Platynereis as a 
powerful lophotrochozoan model for studying neuropep- 
tides. Peptide functions can be easily investigated in 
Platynereis larvae by bath application of synthetic neu- 
ropeptides [28,65]. Recent technical advances now 
allow various genetic and other manipulations of pep- 
tidergic systems. For example, laser nanosurgery [65], 
transgenesis [79], morpholino-mediated gene knock- 
down [65], cellular resolution RNA in situ hybridization 
[80], complete neuronal reconstruction from TEM 
sections [65,81], and whole-body gene expression pat- 
tern registration [66,82] can be applied to explore the 
function of peptidergic systems in Platynereis. 

Conclusion 

Here we used an integrative approach, combining tran- 
scriptomics, peptidomics and sequence homology searches, 
to obtain a broad pNP complement for the model annelid 
Platynereis dumerilii. Using homology-based sequence 
clustering and a comprehensive set of multiple sequence 
alignments and peptide-motif screens, we classified many 
Platynereis pNPs into conserved families. This work will 
serve as the foundation for further research of neuropeptide 



functions in Platynereis and for the study of conserved 
bilaterian peptidergic cells. Our pNP catalog will also pro- 
vide a useful resource for the identification of pNPs in other 
annelids and mollusks (comprising more than 100,000 ex- 
tant species), as well as understudied non-model marine 
invertebrates (e.g. bryozoans or brachiopods). 

Methods 

RNA extraction 

Platynereis samples were obtained from an in-house cul- 
ture at 18°C as previously described [83]. Samples were 
collected for RNA extraction at the following stages: fertil- 
ized egg, 24 hpf, 36 hpf, 48 hpf, 72 hpf, 4 days post 
fertilization (dpf), 10 dpf, 15 dpf, 1 month post fertilization 
with pre-cephalic metamorphosis morphology, 1 month 
fertilization with post-cephalic metamorphosis morph- 
ology, 3 month adult asexual (atokous) worm, male and fe- 
male sexually mature (epitokous) worm. Total larval RNA 
was extracted from pooled batches of larvae (minimum 3 
batches), consisting of hundreds of individuals. Total RNA 
for the 3-month adult stage was extracted from the com- 
bined tissue of 10 worms. Male and female epitoke RNA 
was extracted from the combined tissue of 4 epitokes for 
each sex. 

Total RNA was extracted from all samples using TRIr- 
eagent (Sigma) according to the manufacturer s protocol, 
with two additional phenol:chloroform:isoamyl alcohol 
(25:24:1) phase separations, followed by one additional 
l-Bromo-3-chloropropane (BCP) phase separation step 
prior to isopropanol precipitation. 

Transcriptome sequencing, assembly and annotation 

We used a combination of techniques to obtain a high- 
coverage Platynereis dumerilii transcriptome sequence. 
First, we generated a custom, normalized, full-length, mixed- 
stages cDNA library (with m 7 Gppp affinity purification 
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Figure 8 Clustering of Platynereis pNPs based on stage-specific transcription profiles. Pearson-clustering of Plotynereis pNPs based on 
stage-specific transcription profiles, displayed as normalized reads for each pNP (left) or normalized reads transformed to 1 (right), red: 
high expression, blue: low expression. 
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to avoid bacterial RNA contamination; Invitrogen), cloned 
into the pCMV-Sport6 vector. After plating, we sequenced 
83,152 randomly picked clones using the Sanger technol- 
ogy (ABI 3730) and the SP6 primer. 1,115 of the clones 
were also sequenced with the T7 primer. For base calling 
and vector trimming we used Phred [84,85] and Cross- 
match. Mitochondrial and ribosomal RNA sequences were 
removed using Ribopicker 0.4.3 [86] and a database of 
Platynereis mtDNA and rRNAs. The sequences were 
filtered with dustmasker, and every sequence with more 
than 20% low complexity regions was removed. The result- 
ing 77,419 sequences were submitted to dbEST [GenBank: 
JZ391525 - JZ468943]. Second, we performed 454 sequen- 
cing (GS FLX, Roche/454) on the PCR-amplified cDNA 
library, following concatenation and fragmentation. We ob- 
tained 2,757,258 reads following adaptor trimming, quality 
(0.05), and length filtering (50 bp cutoff), using the soft- 
ware package CLC Genomics Workbench 4.5.1. Third, 
we performed Illumina sequencing on total RNA iso- 
lated from different developmental stages. Construction 
of sample-specific cDNA libraries from 5 \xg total RNA, 
and paired-end transcriptome sequencing with an Illu- 
mina Hiseq 2000, was performed by GATC Biotech 
(Konstanz, Germany). Transcriptome Sequencing data 
was analysed using CLC Genomics Workbench 4.5.1 and 
5.5.1 (CLC Bio). The raw read data for each stage-specific 
library was first filtered to remove Illumina adapter se- 
quences, low quality sequences (Quality Limit 0.05) 
and short fragments (less than 30 base pairs). 

All filtered 454 and Illumina reads were assembled 
using CLC Genomics Workbench 4.5.1. The resulting 
contigs and singletons were joined with all EST se- 
quences, and passed to the CAP3 assembler with default 
parameters [87]. Transcript sequences were searched for 
the longest ORF and translated. 

Transcriptome mapping 

The stage-specific filtered libraries were mapped to a 
subset of the assembled transcriptome of Platynereis (in- 
cluding only those sequences that had a BLASTX hit 
with an e-value < le-5 to the SwissProt database, and the 
98 pNPs, a total of 52,631 transcripts) using the RNA-Seq 
Analysis function, with the following mapping parameters: 
paired distance 250 - 350 base pairs, minimum length 
fraction 0.8, minimum similarity fraction 0.9, maximum 
number of mismatches 2. 

The total number of reads mapped to each gene in 
each stage-specific sample was normalized for total li- 
brary size (reads per million (rpm)). For the subsequent 
analysis we focused only on the 98 pNP genes. To view 
the global pattern of neuropeptide expression through- 
out the Platynereis life cycle, we plotted the total sum of 
normalized expression values for all 98 pNPs in each 
stage-specific library. 



The pNP genes were filtered based on their normal- 
ized expression values using the EdgeR RNA-Seq ana- 
lysis package (Bioconductor) [88] in R version 3.0.1 to 
retain only those genes with at least 2 rpm in at least 
one developmental stage. The expression values were 
also further normalized for gene length (reads per kilo- 
base million (RPKM)). The transcriptional profiles of the 
remaining 79 pNPs were clustered using hierarchical 
clustering with a Pearson correlation distance measure. 
The expression values were plotted using both the nor- 
malized expression values and the same values trans- 
formed to a fraction of 1. 

Peptide extraction 

Platynereis larval and juvenile stages from 2 days post 
fertilization until 15 days post fertilization were collected 
in a small net of 100 \im mesh size (approximately 500- 
1000 animals in total). To remove contamination, ani- 
mals were washed several times in sterile-filtered natural 
seawater. Excessive salts were washed off by rinsing ani- 
mals 1-2 seconds in distilled water. Specimens were 
immediately transferred to ice-cold extraction solution 
(methanol : glacial acetic acid : distilled water, 90:9:1). 
Samples were centrifuged for 1 h at 4°C. The super- 
natant that contained the dissolved peptides was evap- 
orated completely. Peptides were dissolved in 100 \A 
double distilled H 2 0 and used for mass-spectrometry. 

Mass spectrometry: peptide sample preparation and 
LC-MS analysis 

Neuropeptide mixtures were either directly desalted 
with Ci 8 StageTips [89] or reduced and alkylated as de- 
scribed before [90]. LC-MS analyses were performed on 
an EasyLC nano-HPLC (Proxeon Biosystems) coupled 
to an LTQ Orbitrap Elite mass spectrometer (Thermo 
Scientific). Separations of the peptide mixtures were 
done as described elsewhere [91] with slight modifica- 
tions: Peptides were eluted with a 87-min segmented 
gradient of 5-33-90% HPLC solvent B (80% ACN in 
0.5% acetic acid). 

The mass spectrometer was operated in the positive 
ion mode. Precursor ions were recorded in the Orbitrap 
mass analyzer at a resolution of 120,000. The target 
value for the Orbitrap was 10 6 charges and the max- 
imum allowed fill time was 100 ms. The 15 most intense 
precursor ions were sequentially fragmented in each 
scan cycle. High-resolution HCD MS/MS spectra were 
acquired with a resolution of 15,000 and a target value 
of 40,000. The normalized collision energy was set to 35, 
activation time to 0.1 ms and the first mass to 120 
thomson. A minimum of 5000 counts were required to 
trigger MS/MS fragmentation and the maximum allowed 
fill time was 150 ms. The isolation window for MS/MS 
fragmentation was set to 2 thomson. Precursor ions were 
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excluded from sequencing for 60 s after MS/MS. In one of 
the measurements MS/MS on singly charged precursor 
ions was enabled. 

MS data processing, bioinformatic analysis and validation 

In total three MS measurements were performed. The 
acquired MS raw files were processed separately using 
the MaxQuant software (vl.2.2.9) [92]. Extracted peak 
lists were submitted to database search using the An- 
dromeda search engine [93] to query a target-decoy [94] 
database consisting of the predicted Platynereis prote- 
ome (51,767 entries), the predicted secretome (11,075 
entries), predicted neuropeptides (347 entries), com- 
monly observed lab contaminants (248 entries), and the 
reversed complements of a those sequences (130,831 en- 
tries). We required full tryptic specificity allowing up to 
three missed cleavages and set the minimal peptide 
length to five amino acids. The initial precursor mass 
tolerance was set to 6 ppm, for fragment ions we used a 
mass tolerance of 20 ppm. For reduced and alkylated 
samples we defined carbamidomethylation of cysteins as 
fixed modification in the database search. 

In order to screen for potential peptide modifications 
we performed a pre-run on one of the measurements by 
defining a fixed set of commonly observed (oxidation of 
methionine, acetylation of the protein N-terminus) and 
expected variable modifications (amidation of the pep- 
tide C-terminus) and iteratively included the following 
modifications in the database search: amidation of gly- 
cine; methylation of lysine and arginine; sulfation of 
tyrosine; acetylation of serine, threonine, alanine, gly- 
cine; pyroglutamic acid. Based on the results of this 
analysis we chose to include the following variable 
modifications in the actual database search: oxidation 
of methionine, sulfation of tyrosine, acetylation of the 
protein N-terminus, amidation of the peptide C-terminus, 
pyroglumatic acid. 

Identified peptide spectrum matches (PSM) were sta- 
tistically scored by MaxQuant software by calculation of 
posterior error probabilities (PEP) for each PSM. We 
considered all PSMs having a PEP below 0.01 for further 
analysis. The complete list of all identified peptide evi- 
dences including the respective spectra can be found in 
Additional file 9. 

Antibody purification, immunohistology and imaging 

A specific rabbit antibody against CEFLGa was raised and 
affinity-purified following a previously described protocol 
[13]. The specificity of the obtained EFLGa antibody was 
also confirmed by pre-absorbing the antibody in 5 mM 
EFLGa for 2 h before applied in immunostainings, which 
completely blocked the antibody signal. The specificity 
of the cross-species reactive rabbit FVa antibody has 
been documented previously [13]. All samples were 



counterstained with a mouse acetylated tubulin antibody 
(Sigma, Saint Louis, USA). Immunostaining on 48 hpf and 
72 hpf Platynereis larvae as well as image acquisition was 
performed as described [13]. 

Additional files 



Additional file 1: Platynereis transcriptome assembly, part 1 in 
GenBank format (txt). The file was compressed with tar -jcvf archive_- 
name.tar.bz2 file_to_compress, use tar -jxvfto uncompress it. Contains 
85,000 sequences. 

Additional file 2: Platynereis transcriptome assembly, part 2 in 
GenBank format (txt). The file was compressed with tar -jcvf archive_- 
name.tar.bz2 file_to_compress, use tar -jxvfto uncompress it. Contains 
85,000 sequences. 

Additional file 3: Platynereis transcriptome assembly, part 3 in 
GenBank format (txt). The file was compressed with tar -jcvf archive_- 
name.tar.bz2 file_to_compress, use tar -jxvfto uncompress it. Contains 
85,000 sequences. 

Additional file 4: Platynereis transcriptome assembly, part 4 in 
GenBank format (txt). The file was compressed with tar -jcvf archive_- 
name.tar.bz2 file_to_compress, use tar -jxvfto uncompress it. Contains 
96,626 sequences. 

Additional file 5: Platynereis predicted proteins in FASTA format 

(txt). The file was compressed with tar -jcvf archive_name.tar.bz2 file_to_- 
compress, use tar -jxvfto uncompress it. Contains 51,767 sequences. 

Additional file 6: Platynereis predicted secreted proteins in FASTA 
format (txt). The SP was removed from the sequences to facilitate the 
matching of MS hits. Contains 1 1,075 sequences. 

Additional file 7: Platynereis pNPs in FASTA format (txt). Contains 
104 sequences. 

Additional file 8: Precursor structures and peptide logos for 
repetitive motifs for all Platynereis pNPs, in portable document 
format. 

Additional file 9: Mass-spectrometry hits obtained by the analysis 
of Platynereis peptide extracts with the corresponding spectra for 
each hit, in portable document format. 

Additional file 10: Comprehensive list of multiple sequence 
alignments of Platynereis pNPs with homologous sequences, in 
portable document format. 

Additional file 11: Normalized lllumina RNA-sequencing read 
counts for 98 Platynereis pNPs in 13 different life-cycle stages, as 
excel spreadsheet. 



Abbreviations 

AKH: adipokinetic hormone; BLAST: Basic Local Alignment Search Tool; 
CCAP: Crustacean cardio active peptide; conoCAP: Cardio active peptide of 
the conus snail; CRF: Corticotropin releasing factor; DH31: Diuretic hormone 
31; DH44: Diuretic hormone 44; dpf: Days post fertilization; ELH: Egg laying 
hormone; EP: Annelid excitatory peptide; EST: Expressed sequence tag; 
GnRH: Gonadotropin releasing hormone; hpf: Hours post fertilization; 
IRP: Insulin related peptide; LC-MS/MS: Liquid chromatography - tandem 
mass spectrometry; L1 1: Abdominal ganglion neuropeptide L1 1 of A californica; 
MIP: Myoinhibitory peptide; MS: Mass spectrometry; mt: Metatroch; 
NKY: Mollusk neuropeptide KY; NMB: Neuromedin B; nmed-U: Neuromedin U; 
NPF: Neuropeptide F; NPFF: Neuropeptide FF; NPS: Neuropeptide S; 
NPY: Neuropeptide Y; ORF: Open reading frame; patl: First paratroch; 
pat2: Second paratroch; PBAN: Pheromone biosynthesis activating 
neuropeptide; pdf: Pigment dispersing factor; PEP: Posterior error probabilities; 
pNP: Proneuropeptide; PSM: Peptide spectrum matches; pt: Prototroch; 
RPCH: Red pigment concentrating hormone; sCAP: Small cardio active peptide; 
SP: Signal peptide. 

Competing interests 

The authors declare that they have no competing interests. 



Conzelmann et al. BMC Genomics 2013, 14:906 
http://www.biomedcentral.com/1471 -21 64/1 4/906 



Page 13 of 15 



Authors' contributions 

MC extracted peptide samples, analyzed the transcriptome and 
mass-spectrometry data, performed immunoassays, participated in the 
design of the study and drafted the manuscript. EAW generated the 
transcriptome resource, performed the stage-specific transcription profiling and 
helped to draft the manuscript. KK performed the processing, bioinformatic 
analysis and validation of the MS-data and helped to draft the manuscript. 
MFW processed and analyzed mass-spectrometry peptide samples. BM 
designed and coordinated the mass-spectrometry analysis and helped to 
draft the manuscript. GJ designed the study, generated the transcriptome 
resource, performed sequence analysis and sequence clustering and 
drafted the manuscript. All authors read and approved the final 
manuscript. 

Acknowledgements 

This work was supported by Max Planck Society Sequencing Grant 
M.IF.A.EN^A/8050 to GJ. The research leading to these results was supported 
by the European Research Council under European Union Seventh 
Framework Program FP7/2007-2013 and European Research Council Grant 
Agreement 260821. 

Author details 

] Max Planck Institute for Developmental Biology, Spemannstrasse 35, 72076, 
Tubingen, Germany. 2 Proteome Center Tubingen, University of Tubingen, 
Auf der Morgenstelle 1 5, 72076, Tubingen, Germany. 

Received: 12 August 2013 Accepted: 17 December 2013 
Published: 20 December 2013 

References 

1. Liu F, Baggerman G, Schoofs L, Wets G: The construction of a bioactive 
peptide database in metazoa. J Proteome Res 2008, 7:41 1 9-4131. 

2. Watanabe H, Fujisawa T, Holstein TA/V: Cnidarians and the evolutionary 
origin of the nervous system. Dev Growth Differ 2009, 51:167-183. 

3. Srivastava M, Simakov O, Chapman J, Fahey B, Gauthier ME, Mitros T, 
Richards GS, Conaco C, Dacre M, Hellsten U, Larroux C, Putnam NH, Stanke 
M, Adamska M, Darling A, Degnan SM, Oakley TH, Plachetzki DC, Zhai Y, 
Adamski M, Calcino A, Cummins SF, Goodstein DM, Harris C, Jackson DJ, 
Leys SP, Shu S, Woodcraft BJ, Vervoort M, Kosik KS, et al: The Amphimedon 
queenslandica genome and the evolution of animal complexity. Nature 
2010, 466:720-726. 

4. Attenborough RM, Hayward DC, Kitahara MV, Miller DJ, Ball EE: A "neural" 
enzyme in nonbilaterian animals and algae: preneural origins for 
peptidylglycine alpha-amidating monooxygenase. Mol Biol Evol 2012, 
29:3095-3109. 

5. Hook V, Funkelstein L, Lu D, Bark S, Wegrzyn J, Hwang SR: Proteases for 
processing proneuropeptides into peptide neurotransmitters and 
hormones. Annu Rev Pharmacol Toxicol 2008, 48:393-423. 

6. Eipper BA, Stoffers DA, Mains RE: The biosynthesis of neuropeptides: 
peptide alpha-amidation. Annu Rev Neurosci 1992, 15:57-85. 

7. Marks N, Stern F, Kastin AJ: Biodegradation of a-MSH and derived 
peptides by rat brain extracts, and by rat and human serum. Brain 
Res Bull 1976, 1:591-593. 

8. In Y, Ono H, Ishida T: Structural studies on C-amidated amino acids and 
peptides: function of amide group in molecular association in crystal 
structures of Val-Gly-NH2, Ser-Phe-NH2, Gly-Tyr-NH2 and Pro-Tyr-NH2 
hydrochloride salts. Chem Pharm Bull 2002, 50:571-577. 

9. Edison AS, Espinoza E, Zachariah C: Conformational ensembles: the role 
of neuropeptide structures in receptor binding. J Neurosci 1999, 
19:6318-6326. 

10. Merkler DJ: C-terminal amidated peptides: production by the in vitro 
enzymatic amidation of glycine-extended peptides and the importance 
of the amide to bioactivity. Enzyme Microb Technol 1994, 16:450-456. 

11. Prigge ST, Mains RE, Eipper BA, Amzel LM: New insights into copper 
monooxygenases and peptide amidation: structure, mechanism and 
function. Cell Mol Life Sci 2000, 57:1 236-1 259. 

12. Wegener C, Gorbashov A: Molecular evolution of neuropeptides in the 
genus Drosophila. Genome Biol 2008, 9:R1 31 . 

13. Conzelmann M, Jekely G: Antibodies against conserved amidated 
neuropeptide epitopes enrich the comparative neurobiology toolbox. 
Evodevo 2012, 3:23. 



14. Liu F, Baggerman G, Schoofs L, Wets G: Uncovering conserved patterns in 
bioactive peptides in Metazoa. Peptides 2006, 27:3137-3153. 

15. Hewes RS: Neuropeptides and Neuropeptide receptors in the Drosophila 
melanogaster genome. Genome Res 2001, 1 1:1 126-1 142. 

16. Cho HJ, Acharjee S, Moon MJ, Oh DY, Vaudry H, Kwon HB, Seong JY: 
Molecular evolution of neuropeptide receptors with regard to 
maintaining high affinity to their authentic ligands. Gen Comp Endocrinol 
2007, 153:98-107. 

17. Kim D-K, Cho EB, Moon MJ, Park S, Hwang J-l, Do Rego J-L, Vaudry H, 
Seong JY: Molecular coevolution of neuropeptides gonadotropin-releasing 
hormone and kisspeptin with their cognate G protein-coupled receptors. 
Front Neurosci 2012, 6:3. 

1 8. Jurenka R, Nusawardani T: The pyrokinin/pheromone biosynthesis-activating 
neuropeptide (PBAN) family of peptides and their receptors in Insecta: 
evolutionary trace indicates potential receptor ligand-binding domains. 
Insect Mol Biol 201 1 , 20:323-334. 

19. Jekely G: Global view of the evolution and diversity of metazoan 
neuropeptide signaling. Proc Natl Acad Sci U S A 201 3, 1 1 0:8702-8707. 

20. Mirabeau O, Joly J-S: Molecular evolution of peptidergic signaling systems 
in bilaterians. Proc Natl Acad Sci U S A 201 3, 1 1 0:E2028-E2037. 

21. Schoofs L, Holman GM, Hayes TK, Nachman RJ, De Loof A: Isolation, 
identification and synthesis of locusta myoinhibiting peptide (LOM-MIP), 
a novel biologically active neuropeptide from Locusta migratoria. Regul 
Pept 1991, 36:111-119. 

22. Meunier JC, Mollereau C, Toll L, Suaudeau C, Moisand C, Alvinerie P, Butour JL, 
Guillemot JC, Ferrara P, Monsarrat B: Isolation and structure of the endogenous 
agonist of opioid receptor-like ORL1 receptor. Nature 1 995, 377:532-535. 

23. Reinscheid RK, Nothacker HP, Bourson A, Ardati A, Henningsen RA, Bunzow 
JR, Grandy DK, Langen H, Monsma FJ, Civelli O: Orphanin FQ: a 
neuropeptide that activates an opioid like G protein-coupled receptor. 
Science 1995, 270:792-794. 

24. Katafuchi T, Kikumoto K, Hamano K, Kangawa K, Matsuo H, Minamino N: 
Calcitonin receptor-stimulating peptide, a new member of the calcitonin 
gene-related peptide family. Its isolation from porcine brain, structure, 
tissue distribution, and biological activity. J Biol Chem 2003, 278:12046-12054. 

25. Veenstra JA: Neurohormones and neuropeptides encoded by the 
genome of Lottia gigantea, with reference to other mollusks and insects. 
Gen Comp Endocrinol 2010, 167:86-103. 

26. Nathoo AN, Moeller RA, Westlund BA, Hart AC: Identification of 
neuropeptide-like protein gene families in Caenorhabditis elegans and 
other species. Proc Natl Acad Sci 2001, 98:14000-14005. 

27. Hinuma S, Shintani Y, Fukusumi S, lijima N, Matsumoto Y, Hosoya M, Fujii R, 
Watanabe T, Kikuchi K, Terao Y, Yano T, Yamamoto T, Kawamata Y, Habata 
Y, Asada M, Kitada C, Kurokawa T, Onda H, Nishimura O, Tanaka M, Ibata Y, 
Fujino M: New neuropeptides containing carboxy-terminal RFamide and 
their receptor in mammals. Nat Cell Biol 2000, 2:703-708. 

28. Conzelmann M, Offenburger S-L, Asadulina A, Keller T, Munch TA, Jekely G: 
Neuropeptides regulate swimming depth of Platynereis larvae. Proc Natl 
Acad SciUSAlOU, 108:E1 174-E1 183. 

29. McVeigh P, Mair GR, Atkinson L, Ladurner P, Zamanian M, Novozhilova E, 
Marks NJ, Day TA, Maule AG: Discovery of multiple neuropeptide families 
in the phylum Platyhelminthes. Int J Parasitol 2009, 39:1243-1252. 

30. Rowe ML, Elphick MR: The neuropeptide transcriptome of a model 
echinoderm, the sea urchin Strongylocentrotus purpuratus. Gen Comp 
Endocrinol 2012, 179:331-344. 

31. Hummon AB, Amare A, Sweedler JV: Discovering new invertebrate 
neuropeptides using mass spectrometry. Mass Spectrom Rev 2006, 25:77-98. 

32. Boonen K, Landuyt B, Baggerman G, Husson SJ, Huybrechts J, Schoofs L: 
Peptidomics: the integrated approach of MS, hyphenated techniques 
and bioinformatics for neuropeptide analysis. J Sep Sci 2008, 31:427-445. 

33. Menschaert G, Vandekerckhove TTM, Baggerman G, Schoofs L, Luyten W, 
Criekinge WV: Peptidomics coming of Age: a review of contributions 
from a bioinformatics angle. J Proteome Res 2010, 9:2051-2061. 

34. Collins JJ, Hou X, Romanova EV, Lambrus BG, Miller CM, Saberi A, Sweedler 
JV, Newmark PA: Genome-wide analyses reveal a role for peptide 
hormones in planarian germline development. Plos Biol 2010, 8:e1 000509. 

35. Li B, Predel R, Neupert S, Hauser F, Tanaka Y, Cazzamali G, Williamson M, 
Arakane Y, Verleyen P, Schoofs L, Schachtner J, Grimmelikhuijzen CJP, 
Park Y: Genomics, transcriptomics, and peptidomics of neuropeptides 
and protein hormones in the red flour beetle Tribolium castaneum. 
Genome Res 2007, 18:113-122. 



Conzelmann et al. BMC Genomics 2013, 14:906 
http://www.biomedcentral.com/1471 -21 64/1 4/906 



36. Hauser F, Neupert S, Williamson M, Predel R, Tanaka Y, Grimmelikhuijzen G: 
Genomics and peptidomics of neuropeptides and protein hormones present 
in the parasitic wasp Nasonia vitripennis. J Proteome Res 2010, 9:5296-5310. 

37. Dircksen H, Neupert S, Predel R, Verleyen P, Huybrechts J, Strauss J, Hauser 
F, Stafflinger E, Schneider M, Pauwels K, Schoofs L, Grimmelikhuijzen CJP: 
Genomics, transcriptomics, and peptidomics of Daphnia pulex 
neuropeptides and protein hormones. J Proteome Res 201 1, 10:4478-4504. 

38. Xie F, London SE, Southey BR, Annangudi SP, Amare A, Rodriguez-Zas SL, 
Clayton DF, Sweedler JV: The zebra finch neuropeptidome: prediction, 
detection and expression. BMC Biol 2010, 8:28. 

39. Oumi T, Ukena K, Matsushima 0, Ikeda T, Fujita T, Minakata H, Nomoto K: 
Annetocin - an oxytocin-related peptide isolated from the earthworm, 
Eisenia foetida. Biochem Biophys Res Commun 1 994, 1 98:393-399. 

40. Veenstra JA: Neuropeptide evolution: neurohormones and neuropeptides 
predicted from the genomes of Capitella teleta and Helobdella robusta. 
Gen Comp Endocrinol 201 1 , 1 71 :1 60-1 75. 

41. Salzet M, Bulet P, Wattez C, Malecha J: FMRFamide-related peptides in the 
sex segmental ganglia of the Pharyngobdellid leech Erpobdella 
octoculata. Identification and involvement in the control of hydric 
balance. Eur J Biochem 1994, 221:269-275. 

42. Evans BD, Pohl J, Kartsonis NA, Calabrese RL: Identification of RFamide 
neuropeptides in the medicinal leech. Peptides 1991, 12:897-908. 

43. Simon TW, Schmidt J, Calabrese RL: Modulation of high-threshold transmission 
between heart interneurons of the medicinal leech by FMRF-NH2. 

J Neurophysiol 1 994, 71 :454-466. 

44. Diaz-Miranda L, de Motta GE, Garcia-Arraras JE: Monoamines and 
neuropeptides as transmitters in the sedentary polychaete Sabellastarte 
magnified: Actions on the longitudinal muscle of the body wall. J Exp 
Zool 1992, 263:54-67. 

45. Krajniak KG, Greenberg MJ: The localization of FMRFamide in the nervous 
and somatic tissues of Nereis virens and its effects upon the isolated 
esophagus. Comp Biochem Physiol 1992, 101:93-100. 

46. Ukena K, Oumi T, Matsushima O, Takahashi T, Muneoka Y, Fujita T, Minakata H, 
Nomoto K: Inhibitory pentapeptides isolated from the gut of the earthworm, 
Eisenia foetida. Comp Biochem Physiol A Physiol 1996, 1 14:245-249. 

47. O'Gara BA, Brown PL, Dlugosch D, Kandiel A, Ku JW, Geier JK, Henggeler NC, 
Abbasi A, Kounalakis N: Regulation of pharyngeal motility by FMRFamide 
and related peptides in the medicinal leech, Hirudo medicinalis. Invert 
Neurosci 1999, 4:0041-0053. 

48. Herbert Z, Pollak E, Zougman A, Boros A, Kapan N, Molnar L: Identification 
of novel neuropeptides in the ventral nerve cord ganglia and their 
targets in an annelid worm, Eisenia fetida. J Comp Neurol 2009, 514:415-432. 

49. Matsushima O, Takahashi T, Morishita F, Fujimoto M, Ikeda T, Kubota I, Nose T, 
Miki W: Two S-lamide Peptides, AKSGFVRIamide and VSSFVRIamide, Isolated 
from an Annelid, Perinereis vancaurica. Biol Bull 2002, 184:216-222. 

50. Jekely G, Colombelli J, Hausen H, Guy K, Stelzer E, Nedelec F, Arendt D: 
Mechanism of phototaxis in marine zooplankton. Nature 2008, 456:395-399. 

51 . Minakata H, Fujita T, Kawano T, Nagahama T, Oumi T, Ukena K, Matsushima O, 
Muneoka Y, Nomoto K: The leech excitatory peptide, a member of the 
GGNG peptide family: isolation and comparison with the earthworm GGNG 
peptides. FEBS Lett 1997, 410:437-442. 

52. Matsushima O, Takahama H, Ono Y, Nagahama T, Morishita F, Furukawa Y, 
Iwakoshi-Ukena E, Hisada M, Takuwa-Kuroda K, Minakata H: A novel 
GGNG-related neuropeptide from the polychaete Perinereis vancaurica. 
Peptides 2002, 23:1379-1390. 

53. Nagahama T, Ukena K, Oumi T, Morishita F, Furukawa Y, Matsushima O, 
Satake H, Takuwa K, Kawano T, Minakata H, Nomoto K: Localization of leech 
excitatory peptide, a member of the GGNG peptides, in the central 
nervous system of a leech {Whitmania pigra) by immunohistochemistry 
and in situ hybridization. Cell Tissue Res 1999, 297:155-162. 

54. Salzet M, Verger-Bocquet M, Vandenbulcke F, Van Minnen J: Leech egg- 
laying-like hormone: structure, neuronal distribution and phylogeny. 
Brain Res Mol Brain Res 1 997, 49:21 1 -221 . 

55. Wang Y, Price DA, Sahley CL: Identification and characterization of a 
myomodulin-like peptide in the leech. Peptides 1998, 19:487-493. 

56. Tobin A-E, Calabrese RL: Myomodulin increases Ih and inhibits the Na/K 
pump to modulate bursting in leech heart interneurons. J Neurophysiol 
2005, 94:3938-3950. 

57. Takahashi T, Matsushima O, Morishita F: A myomodulin-CARP-related 
peptide isolated from a polychaete annelid, Perinereis vancaurica. Zoolog 
Sci 1994, 11:33-38. 



Page 14 of 15 



58. Tessmar-Raible K, Raible F, Christodoulou F, Guy K, Rembold M, Hausen H, 
Arendt D: Conserved sensory-neurosecretory cell types in annelid and 
fish forebrain: insights into hypothalamus evolution. Cell 2007, 

129:1389-1400. 

59. Wagenaar DA, Hamilton MS, Huang T, Kristan WB, French KA: A 
hormone-activated central pattern generator for courtship. Curr Biol 
2010, 20:487-495. 

60. Tsai PS, Zhang L: The emergence and loss of gonadotropin-releasing 
hormone in protostomes: orthology, phylogeny, structure, and function. 

Biol Reprod 2008, 79:798-805. 

61 . Luis Quintanar J, Gutierrez-Garcia K, Castillo-Hernandez L: The neuropeptide 
Gonadotropin-releasing hormone modifies the spontaneous muscular 
contraction in the earthworm: Eisenia foetida. Invert Neurosci 201 1 , 11:11 3-1 1 6. 

62. LeRoith D, Lesniak MA, Roth J: Insulin in Insects and Annelids. Diabetes 
1981,30:70-76. 

63. Diaz-Miranda L, Escalona De Motta G, Garc A-A, JE S: Localization of neuropeptides 
in the nervous system of the marine annelid Sabellastarte magnifica. Cell Tissue 
Res 1991,266:209-217. 

64. Curry WJ, Fairweather I, Johnston CF, Halton DW, Buchanan KD: 
Immunocytochemical demonstration of vertebrate neuropeptides in the 
earthworm Lumbricus terrestris (Annelida, Oligochaeta). Cell Tissue Res 
1989, 257:577-586. 

65. Conzelmann M, Williams EA, Tunaru S, Randel N, Shahidi R, Asadulina A, 
Berger J, Offermanns S, Jekely G: Conserved MIP receptor-ligand pair 
regulates Platynereis larval settlement. Proc Natl Acad Sci U S A 201 3, 
110:8224-8229. 

66. Tomer R, Denes AS, Tessmar-Raible K, Arendt D: Profiling by image 
registration reveals common origin of annelid mushroom bodies and 
vertebrate pallium. Cell 2010, 142:800-809. 

67. Fischer AH, Hen rich T, Arendt D: The normal development of Platynereis 
dumerilii (Nereididae, Annelida). Front Zool 2010, 7:31. 

68. Fischer A, Dorresteijn A: The polychaete Platynereis dumerilii (Annelida): a 
laboratory animal with spiralian cleavage, lifelong segment proliferation 
and a mixed benthic/pelagic life cycle. Bioessays 2004, 26:314-325. 

69. Petersen TN, Brunak S, Heijne Von G, Nielsen H: SignalP 4.0: discriminating 
signal peptides from transmembrane regions. Nat Methods 201 1, 8:785-786. 

70. Mbikay M, Seidah NG, Chretien M: Neuroendocrine secretory protein 7B2: 
structure, expression and functions. Biochem J 2001, 357:329-342. 

71. Senatorov W, Yang CR, Marcinkiewicz M, Chretien M, Renaud LP: 
Depolarizing action of secretory granule protein 7B2 on Rat supraoptic 
neurosecretory neurons. J Neuroendocrinal 1993, 5:533-536. 

72. Rayne RC, O'Shea M: Reconstitution of adipokinetic hormone biosynthesis 
in vitro indicates steps in prohormone processing. Eur J Biochem 1994, 
219:781-789. 

73. Frickey T, Lupas A: CLANS: a Java application for visualizing protein 
families based on pairwise similarity. Bioinformatics 2004, 20:3702-3704. 

74. Baggerman G, Liu F, Wets G, Schoofs L: Bioinformatic analysis of peptide 
precursor proteins. Ann N Y Acad Sci 2005, 1040:59-65. 

75. Aravind L, Koonin EV: A colipase fold in the carboxy-terminal domain of 
the Wnt antagonists - the Dickkopfs. Curr Biol 1998, 8:R477-R478. 

76. Dunn JG, Foo CK, Belletier NG, Gavis ER, Weissman JS: Ribosome profiling 
reveals pervasive and regulated stop codon readthrough in Drosophila 
melanogaster. eLife 2013, 2:e01 179. 

77. Tessmar-Raible K, Arendt D: Emerging systems: between vertebrates and 
arthropods, the Lophotrochozoa. Curr Opin Genet Dev 2003, 13:331-340. 

78. Robb SMC, Ross E, Sanchez Alvarado A: SmedGD: the Schmidtea 
mediterranea genome database. Nucleic Acids Res 2008, 36:D599-D606. 

79. Backfisch B, Veedin Rajan VB, Fischer RM, Lohs C, Arboleda E, Tessmar-Raible 
K, Raible F: Stable transgenesis in the marine annelid Platynereis dumerilii 
sheds new light on photoreceptor evolution. Proc Natl Acad Sci USA 
2013, 110:193-198. 

80. Jekely G, Arendt D: Cellular resolution expression profiling using confocal 
detection of NBT/BCIP precipitate by reflection microscopy. Biotech 2007, 
42:751-755. 

81. Randel N, Bezares-Calderon LA, Guhmann M, Shahidi R, Jekely G: Expression 
dynamics and protein localization of rhabdomeric opsins in Platynereis 
larvae. Integr Comp Biol 2013, 53:7-16. 

82. Asadulina A, Panzera A, Veraszto C, Liebig C, Jekely G: Whole-body gene 
expression pattern registration in Platynereis larvae. Evodevo 2012, 3:27. 

83. Hauenschild CFA: Platynereis dumerilii. Mikroskopische Anatomie, 
Fortpflanzung, Entwicklung. Stuttgart: Gustav Fischer; 1969. 



Conzelmann et al. BMC Genomics 2013, 14:906 
http://www.biomedcentral.com/1471 -21 64/1 4/906 



Page 15 of 15 



84. Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated 
sequencer traces using phred. I. Accuracy assessment. Genome Res 1998, 
8:175-185. 

85. Ewing B, Green P: Base-calling of automated sequencer traces using 
phred. II. Error probabilities. Genome Res 1998, 8:186-194. 

86. Schmieder R, Lim YW, Edwards R: Identification and removal of ribosomal 
RNA sequences from metatranscriptomes. Bioinformatics 2012, 28:433-435. 

87. Huang X, Madan A: CAP3: A DNA sequence assembly program. Genome 
Res 1999, 9:868-877. 

88. Robinson MD, McCarthy DJ, Smyth GK: edgeR: a Bioconductor package for 
differential expression analysis of digital gene expression data. 
Bioinformatics 2010, 26:139-140. 

89. Rappsilber J, Mann M, Ishihama Y: Protocol for micro-purification, 
enrichment, pre-fractionation and storage of peptides for proteomics 
using StageTips. Nat Protoc 2007, 2:1896-1906. 

90. Borchert N, Dieterich C, Krug K, Schutz W, Jung S, Nordheim A, Sommer RJ, 
Macek B: Proteogenomics of Pristionchus pocificus reveals distinct 
proteome structure of nematode models. Genome Res 2010, 20:837-846. 

91. Franz-Wachtel M, Eisler SA, Krug K, Wahl S, Carpy A, Nordheim A, 
Pfizenmaier K, Hausser A, Macek B: Global detection of protein kinase 
D-dependent phosphorylation events in nocodazole-treated human 
cells. Mol Cell Proteomics 2012, 1 1:160-1 70. 

92. Cox J, Mann M: MaxQuant enables high peptide identification rates, 
individualized p.p.b.-range mass accuracies and proteome-wide protein 
quantification. Nat Biotechnol 2008, 26:1367-1372. 

93. Cox J, Neuhauser N, Michalski A, Scheltema RA, Olsen JV, Mann M: 
Andromeda: a peptide search engine integrated into the MaxQuant 
environment. J Proteome Res 201 1, 10:1794-1805. 

94. Elias JE, Gygi SP: Target-decoy search strategy for increased confidence in 
large-scale protein identifications by mass spectrometry. Nat Methods 
2007, 4:207-214. 



doi:1 0.1 186/1471-2164-14-906 

Cite this article as: Conzelmann et al:. The neuropeptide complement of 
the marine annelid Platynereis dumerilii. BMC Genomics 2013 14:906. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at f~\ RiftMM i rpntral 

www.biomedcentral.com/submit momea central 



